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ABSTRACT 

By investigating the canonical commutation rules for gravitating quantized particles 
in a 2+1 dimensional world it is found that these particles live on a space-time lattice. The 
space-time lattice points can be characterized by three integers. Various representations 
are possible, the details depending on the topology chosen for energy-momentum space. 
We find that an S2 x Si topology yields a physically most interesting lattice within which 
first quantization of Dirac particles is possible. An S3 topology also gives a lattice, but 
does not allow first quantized particles. 



1. INTROUCTION 



Reducing the number of space dimensions from three to two in quantum gravity is a 
quite severe simplification that drastically changes the physics one attempts to describe. 
Due to the absence of local dynamical degrees of freedom Einstein's equations for gravity 
become almost trivial: in the absence of matter space-time is featureless (flat, if there is 
no cosmological constant). Hence pure gravity is a topological theory in this case, and it 
has been investigated by a number of authors 1 . 

Pure 2+1 dimensional gravity is of interest if one wishes to contemplate cosmologies 
with highly non-trivial boundary conditions in the space-time coordinates. However if 
one wishes to gain insights in local aspects of gravity these boundary conditions are not 
so important. The theory then becomes non-trivial only if one adds matter. One of the 
simplest forms of matter a particle physicist can imagine is a single, non-interacting scalar 
field with the following Lagrangian: 

^matter = ^""^^ " W ^ ) • ( L1 ) 

In the absense of gravity the quantum version of this theory generates a Hilbert space 
spanned by all N particle states, where each particle has a given momentum, and the 
particles do not interact. 

What gravity does to classical point particles in 2+1 dimensions is also clear : each 
particle cuts a cone in 2-space, as seen from its own Lorentz rest frame. Since the sur- 
rounding space-time is completely flat this system is still exactly solvable, although in the 
iV particle case the book-keeping of the transitions after a long stretch of time can become 
infinitely complicated and chaotic phenomena can easily occur 3 . 

Does this suggest that the quantum mechanical version of the gravitating point par- 
ticle theory should also be soluble? Many classically integrable systems are also quantum 
mechanically integrable. But it appears that we can compare the classical theory best with 
that of classical particles in a triangular box with arbitrary angles, which is classically inte- 
grable only for finite time intervals, but at large times chaotic, and quantum mechanically 
not exactly solvable. 

Actually, the situation for our toy model is much worse than for the particle in a 
triangular box: we do not really know how to define the corresponding quantum model. 
Is it a functional integral? If we try to define the theory from that angle we must face the 
fact that it is not renormalizable, so that one must seriously doubt whether the functional 
integral can be at all meaningfully defined. Should we build a Hilbert space from states 
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with N (moving) punctures in space and time? But that is a first quantized theory. If 
we want a theory that in the absense of matter would approach a scalar field theory we 
want something that is second quantized. This would require the introduction of creation 
and annihilation operators. Creating or annihilating punctures is tricky if we want to use 
coordinates to localize these operators, since they backfire on the coordinates themselves. 

A promising avenue may be first to construct a Hilbert space containing N Dirac like 
gravitating particles, obeying Fermi-Dirac statistics or at least some sort of exclusion 
principle. This would be the first quantized model. Then one might hope to be able 
to construct states where the Dirac sea of negative-energy states is filled, and introduce 
second quantization that way. This programme has not yet been carried out completely. 
A problem is that the negative-energy particles in the Dirac sea would generate negative 
curvature which should somehow have to be compensated, perhaps by a cosmological 
constant. 

In spite of all these obstacles it is of importance to try to find our way in this jungle. 
If a meaningful model can be constructed it might show us on the one hand how gravity 
takes care of its own renormalization problems in the ultraviolet region, and on the other 
hand, in the infrared region, how to handle the system as a quantum cosmology, which 
requires a definition of cosmological time, and a way to handle the measurement problem. 

It has been noted before that the time coordinate must probably be chosen to be a dis- 
crete multiple of a Planckian time unit 5 . This important observation follows from the fact 
that the Hamiltonian is an angular variable, well defined only modulo 2tt (in conveniently 
chosen Planck units). Time is also directly linked to angular momentum. But what about 
the spacelike coordinates? This author's first attempts to identify appropriate variables 
for the iV-particle Hilbert space did not lead to space quantization . This was because 
we used the so-called polygon representation 3 ' 5 ' 6 ' 7 . Natural parameters to characterize an 
iV-particle state appeared to be the lengths Lj of the edges of polygons used to tesselate 
a two-dimensional Cauchy surface. The angles of the polygons are fixed as soon as the 
Lorentz boost parameters are given; that is, the Lorentz boosts rji relating the coordinate 
frame of one polygon to that of a neighboring polygon to which it is attached. At each Li 
we have one rji, and it was discovered that the Lj and r\i are conjugated variables: 



{2^,L,} = Sij. (1.2) 

But the "momenta" rji are not angles, they are hyperbolic angles. The equations 
of motion relate the trigonometric functions of the Hamiltonian, sinH and cosi/, to the 
hyperbolic sines and cosines of the rj parameters. If anything, this would suggest that the 
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imaginary parts, not the real parts, of the parameters Li are to be chosen discrete. If the 
Li were to represent the space coordinates, they still appeared to live in a continuum. 

But using the poygon representation to formulate Hilbert space leads to other problems 
as well. We have very complicated boundary conditions in the form of monodromies as 
soon as particles make a full swing around each other. Also there are constraints: each 
polygon must be a closed polygon, which appears to correspond to the constraint that 
the state must be invariant under Lorentz transformations of the coordinate frames inside 
each polygon. The fact that the exterior angles of each polygon must add up to 2n 
corresponds to the constraint of invariance with respect to time boosts for the frame inside 
each polygon. We found that these constraints are easier to deal with if instead of the Li 
we use coordinates {xi, yi} to localize each corner of each polygon in its own frame. This 
required to deduce the equations of motion and the Poisson brackets from scratch once 
again. 

It is here that we hit upon a pleasant surprise: the momenta conjugated to these Xi 
and yi are true angles once again, as if space itself were limited to discrete numbers as 
well. Is space as well as time a lattice? If this is so it will be hard to recover rotational 
and Lorentz invariance. Our first attempts were to replace the wave equations by lattice 
equations which reduce Poincare invariance to one of its discrete subgroups. But this 
seriously jeopardized the entire structure of the theory. If our monodromies were to be 
restricted to lattice group elements it would become impossible to consider continuum 
limits. This wasn't the theory we were looking for. 

The present paper explains how to set up a more satisfactory scheme. Complete 
invariance under the continuum Lorentz group needs not be sacrificed even if space and 
time are restricted to points on a lattice. Everything is derived from the Poisson algebra, 
which leads one to conclude that the results are inevitable. It is somewhat puzzling to 
the present author why the lattice structure of space and time had escaped attention from 
other investigators up till now. 

There are different possible representations of the lattice. We first thought that the S3 
topology for energy-momentum space would be the most natural choice. It will turn out 
however that in the lattice thus obtained no Lorentz-invariant quantum particles can exist 
unless negative probabilities are accepted. Even a Dirac square root of the partial difference 
equations for the wave functions does not allow for a positive probability interpretation. 
We will demonstrate that this disease can be removed if we use an £2 x Si topology instead. 
The lattice looks only slightly different, and the Dirac equation becomes a beautiful first 
order difference equation. 
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2. THE EQUATIONS FOR ONE PARTICLE 



For a particle at rest, surrounding space is a cone with a deficit angle (3 that can 
be identified with the mass m. When the particle moves it is convenient to consider the 
relation between half the deficit angles at motion and at rest. Therefore we will use units 
such that the deficit angle, which is easily seen to be additively conserved, is identified 
with twice the Hamiltonian H of the system. We will use the symbol \x for half the deficit 
angle at rest, being a parameter that corresponds to mass of a particle. For the time being 
we will assume that \i < n/2. Let the particle move with velocity v = tanh£, where £ is 
the particle's Lorentz boost parameter. The geometry is then sketched in Figure 1. 




Fig. 1. Wedge cut out by moving particle (dot). £ is the boost parameter for the veloc- 
ity of the particle, r\ the one for the velocity of one edge of the wedge. The hamiltonian 
H is half the wedge angle. 

By choosing the wedge such that the velocity vector lies in the direction of its bisectrix 
one achieves that the jump across the wedge is purely spacelike, and gluing the edges 
together provides with a Cauchy surface. Lorentz contraction causes the wedge of a particle 
to widen when it moves. The velocity of the edges, tanh/y, easily follows from the geometry. 
The result is 2 ' 3 ' 5 ' 6 : 



tanH = cosh£ tan^; (2.1) 

tanh/y = sin if tanh£; (2.2) 

cos \i = cos H cosh rj ; (2.3) 

sinh?7 = sin// sinh£. (2-4) 

As explained in the introduction, since r\ is canonically conjugated to the length L of the 
wedge, the hyperbolic sines and cosines suggest a discretization of the imaginary part of 
L, a somewhat mysterious feature. 

When more particles are being considered the total Hamiltonian will be the sum of 
the individual contributions, but there will also be "gravitational" contributions from the 
vertex points between triples of polygons. This is because the Cauchy surface is curved 
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at these points, so that the angles there do not add up to 2-7T. These contributions are 
completely determined by the 77's, as in Eq. (2.3), but we will not discuss them further in 
this paper. 

Now consider the coordinates (x,y) of the particle. First let's have the particle move 
in the ^-direction. Let p x be the canonical momentum associated to x. What is p x l One 
must require 

^-x = tanh£ = |5. (2.5) 
dt dp x 

This is a partial derivative where \i is kept fixed. The Hamiltonian has been postulated to 
be half the wedge angle. We use Eq. (2.3) to derive 

d(cos H cosh 77) = — > 

tanHdH = tanh??d?7 — > 
dH sinHdH 



dp. 



tanh £ tanh rj 

COS [1 

cosHdn = — - — drj = (cos /j)d arctan(sinh?y) 
cosh rj 



provided we keep 



(2.6) 



where also Eq. (2.2) was used. 
This we can write as 

p x = 6 cosn > tan6> = sinh?], (2.7) 

which is synonymous to 

sin6> = tanhry ; (2.8) 

cos6> cosh?? = 1 ; (2.9) 

tan(±0) = tanh(i?7) , (2.10) 



|0| < (2.11) 



(a condition that we will ignore later). Thus we find that, apart from a factor cos/U, the 
canonical momentum p x is an angle 9. 
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3. ENERGY-MOMENTUM SPACE, AND THE SPACE-TIME LATTICE 



Our first impression was that therefore the coordinate x must be quantized, 



X = TL\j COS fl , 



(3.1) 



and similarly the coordinate y. But then one hits two problems: first, space would become 
a rectangular lattice, so that invariance under continuous rotations would be lost, and 
furthermore, if the particle does not move in the x-direction, Eqs. (2.6) — (2.11) would 
hold only for the absolute value \p\, not for p x and p y separately. We have to look at 
momentum space more carefully. There are two angles: the value 9 = pj cos/U, and the 
angle ip describing the direction in which the particle moves, which must be the direction 
of the vector p. This, of course, is the property of line segments lying on an S2 sphere. 
Let us therefore write 



Then (Qi,Q2,Qs) is a three-vector whose length is immaterial, whose angle with the 3- 
axis corresponds to the total momentum, and whose angle in the 1-2 plane is the angle of 
motion of the particle. 

The space coordinates x and y should canonically conjugated^ to these variables p x 
and p y . Since p x and p y together form the points of an S2 sphere it is obvious what 
these conjugated variables are: they are the quantum numbers (£, m) of the spherical har- 
monics. These indeed form a two-dimensional lattice, but one not quite as simple as a 
rectangular one. Note that they allow for the full rotation group U(l) as well as discrete 
translations - the latter being obtained by multiplying a spherical harmonic with another 
spherical harmonic, using the Clebsch-Gordan coefficients for the addition of angular mo- 
menta. How does this lattice have to be combined with the timelike lattice, and how do 
we perform Lorentz transformations? 

The relation between the Hamiltonian and the momentum (that is, the Schrodinger 
equation), can be read off from Eqs. (2.9) and (2.3): 



tan^ / cos fi) =Qi/Q^; 
tan(p < y I cos fi) =Q2/Qd,- 



(3.2) 



cos H = cos fi cos 6 . 



(3.3) 



t Note however the remark at the end of Sect. 3. 
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Furthermore we have 



tan 9 = sin \x sinh £ ; (3.4) 
sin 9 = sin H tanh £ . (3.5) 

These equations are to be compared with the usual relations for non-gravitating relativistic 
particles (corresponding to the ji — » limit): 



H 2 =l x 2 +p 2 - (3.3a) 

[XV 

p = Hv . (3.5a) 



p = -if= ; (3.4*) 



Let us now investigate how a Lorentz transformation in the x direction affects p x and 
H. Let us increase the boost parameter £ by a small Lorentz boost e: 



We find from Eq. (2.7): 



and, using Eq. (2.1), 



Clearly, the couple 



sinh £ — > sinh £ + £ cosh £ . 



tan — > tan 6* + £ sin ^ cosh £ 

= tan 6> + e(cos fi) tan ff , 



tan H — > tan H + e tan |it sinh £ 

cos 



tan6> 
cos ji tan iif / ' 



which we can also write as 



(3.6) 



(3.7) 



e n (3-8) 

tan#H tan6>. v ' 



(3.9) 
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sin fi 



( 



sinh£ 
cosh£ 



) 



(3.9a) 



transforms as a Lorentz vector. One may write 




(3.10) 



where Q4 is again a real quantity, so that we have invariance under the Lorentz transfor- 
mation 



where L is an ordinary 50(2, 1) matrix. Q 3 is not involved in Lorentz transformations. 
There now are topologically distinct ways to represent energy-momentum space. The 
length of the Q vector can be chosen freely. Since we are dealing with angles, not hyperbolic 
angles, it is not advised to use hyperbolic spaces but rather compact spaces. If we take 
Ql 2 + Q\ + Q\ + Q\ = 1, this space forms a Euclidean £3 sphere. The SO (4) spherical 
harmonics then generate a lattice in space and time. Our first investigations involved this 
lattice, and the field equations on it. We now defer further discussion of this lattice to the 
Appendix, because it was discovered that the field equations one then obtains are difference 
equations but they link five consecutive time layers, which prohibits the construction of any 
quantity that can serve as a conserved, non-negative probability density. The square roots 
of these equations, a discrete analogy of the Dirac equation, still links three consecutive 
time layers, which is still physically not acceptable. We explain this situation further in 
the Appendix. 

In our search for a physically most appealing theory we found that the £3 topology 
for energy- momentum space must be abandoned, replacing it by the S2 x Si torus. Thus 
we keep the £2 sphere for momentum space and put the Hamiltonian on Si. In this 
energy-momentum space Lorentz transformations are still possible. We keep the 3-vector 



(where cp is the angle in which the particle moves) as being the one that transforms as 
usual under Lorentz transformations, but now define points on the space-time lattice to 




(3.11) 




(3.12) 
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correspond to harmonic functions on this torus: 

(£,m,t) <- Y im (e^)e- itH , (3.13) 

where t is an integer denoting time. It is clear that here £ roughly corresponds to the 
distance from the origin*, whereas m exactly corresponds to angular momentum in the 
2-plane. A discrete Fourier transformation provides something like an angular position 
with respect to the origin, so that the picture that emerges is as given in Fig. 2a. Depicted 
there is the case when £ takes integral values. Of course we also expect the possibility of 
half-integral values, in which case the space lattice is as sketched in Fig. 2b. Since then 
angular momentum has half-odd-integral values we expect this case to be appropriate for 
fermions. We have not contemplated the possibility of having 'anyons' in this theory (see 
some remarks in the Discussion Section). 




Fig. 2. Lattice in 2-space, a) the bosonic case, b) for fermions. 

It should be emphasized that a Lorentz transformation will not correspond to a simple 
point transformation on our lattice. Lorentz transformations do transform the harmonic 
functions into linear superpositions of such functions. This implies that a Lorentz transfor- 
mation is not a local transformation in this theory. It is of course a point transformation 
in energy-momentum space. 

A second remark is that the coordinates x and y most closely correspond to the angular 
momentum operators L2 and L\ on the 2-sphere: 

£ = — ; 2/=—^; L = L 3 . (3.14) 

COS jJL COS |U 

* More precisely: £{£ + 1) = (x 2 + y 2 ) / cos 2 ^ + m 2 . 
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Consequently, we have the commutation rules 



L ^' yj "cos 2 ^' 

[L,x]=iy, (3 ' 15) 
[L,y] = -ix. 

Here we reinserted the factors cos/U from Eq. (2.7). 

A third remark: the commutation rules between x and p appear to deviate from the 
standard expressions. This is because the momenta have become angles. The commutation 
rules are now the ones generated by the spherical harmonics Yg m in the usual way. This 
is a delicate feature since we used Eq. (2.5) as a starting point. It may be verified that 
Eq. (2.5) holds in the limit of large distances and large times, which is what one really 
must require for the classical limit. 

4. A FIELD EQUATION 

Let us consider the 2-sphere 

Ql + Ql + Ql = i, (4.1) 

and define the momenta as in Eq. (3.2). The mass shell condition, Eq. (3.3), can then be 
rewritten as 

e iH + e -m = (2oos//)<9 3 . (4.2) 

It is not difficult to see what this equation implies for functions defined on our space-time 
lattice. First consider any function ip(£,m,t). Since the functions H and Qi are defined 
on its 'Fourier transform' 

^ (0, <p,H)= Y im(0, <p)e iHt m m, t) , (4.3) 

£,m,t 

the functions e ±lH and Qi act as operators on ip. The first cause a shift in time by one 
step down and up, respectively. The action of the operators Qi is found as follows: 



11 



(Qi + 1Q2) \£, m) = C{£) y / (£ + m + l)(£ + m + 2) \£ + 1, m + 1) 
- C{£ - l)y/(£ - m){£ - m - 1) \£ - 1, m + 1) , 

Qsl^m) = C{£) V(£ + m + 1)(£ - m + 1) \£ + 1, m) 
+ C(£ - 1) v 7 ^ + m)(£ - m) \£ - 1, m), 

(Qi-*Q 2 )K, m> = -C(£) - m + 1)(£ -m + 2)\£ + l,m-l) 
+ C{£ - 1) v 7 ^ + m){£ + m-l)\£-l, m - 1) . 

where C(£) is a normalization factor, depending only on £. These equations are derived us- 
ing standard Clebsch-Gordan techniques: one considers the Q functions in these equations 
as being proportional to the states |1, 1), |1,0) and |1, —1), respectively. The coefficients 
are then found by considering the operators L± on these states. The normalization factor 
C(£) is found, after a series of algebraic manipulations, by inserting Eq. (4.1): 

C{£) = ((2£ + l)(2£ + 3))" 1/2 . (4.7) 

Inserting Eq. (4.5) into Eq. (4.2) now yields a difference equation which in the con- 
tinuum limit should reproduce the Klein-Gordon equation. Notice that this difference 
equation only involves nearest neighbors in the {£, m, t} lattice: 



(4.4) 
(4.5) 
(4.6) 



ip(£, m,t+l) + ip(£, m, t - 1) = (2 cos 



+ m+l)(£-m+l) 
(2£+l)(2£ + 3) 



*/>(£+!, m,t) 



+ 



(£ + m)(£-m) 



(2£- 1)(2£ + 1) 
For large £, small m and small \i this approaches to: 



ip(£-l, m,t) 



(4.8) 



ip(£, m,t+l)- 2ip(£, m, t) + if>(£, m,t-l) = 



ip(£ + 1, m, t) - 2ip(£, m, t) + ip(£ - 1, m, t) 



2- ^ 



m 2 1 
~ W + 8£ 2 i 



x 



2 1 

2 , m ~ 4 

n + 



£ 2 



tp(£, m, t) . 



(4.9) 



Identifying here m with angular momentum and £ with r 
Klein-Gordon equation in polar coordinates. 



\J x 2 + y 2 , we recognise the 
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The fact that Eq. (4.8) emerges as a difference equation involving only nearest neigh- 
bors is the primary physical reason for introducing the {£, m, t} lattice as proposed in this 
paper. 

5. THE LATTICE DIRAC EQUATION 

Since Eq. (4.8) connects three time slices it is not suitable for serving as a quantum 
wave function. We believe that a promising route towards a useful quantized model for 
gravitating particles in 2+1 dimensions is to follow Dirac's philosophy: find an acceptable 
wave function for fermions, which allows for an interpretation of \ip\ 2 as a probability 
distribution, after which we handle the negative energy solutions by filling all these levels. 
This paper will not go that far (there are various important problems that would yet have 
to be addressed), but the first step, a Lorentz covariant equation with a non negative 
preserved probability distribution, can now be made. 

Instead of Eq. (4.2) we need an equation containing e~ lH only, so that an evolution 
operator would be obtained connecting exactly two time slices. To do this, we need not 
only an expression for cos if, as in Eq. (3.3), but also one for sinH. We have, on the sphere 
(4.1): 

sin 2 H = sin 2 ^ + cos 2 \x sin 2 6> = 
= sin 2 /x + cos 2 // (Ql + Ql) 

Therefore 



e lH = cos/uQs — n/sin 2 ^ + cos 2 ^ {Q\ + Ql) = 

I (5-2) 

= cos [iQs - iyQl + Ql + sin 2 ^Q| . 

To have the equation homogeneous in the Q's will be of importance later. Clearly, the 
square root in this equation would turn it into a non-local one on the lattice. Following 
Dirac we replace this equation by a linear one: 

e~ tH tfj = [cos /j,Q s -iaiQi -ia^Qi - i/? sin ^ Q 3 ] ^ , (5.3) 

where and (3 obey the usual Dirac anticommutation rules, so that the Eigenvalues of 
the operator in (5.3) correspond to the field equation (4.8). 

13 



(5-1) 



To check Lorentz invariance we consider the Lorentz vector (3.11), (3.12). 
Q4, as defined in (3.10), obeys 



Ql + Ql = 1/cosV (5.4) 

so that 

e~ iH = cos/i(Q 3 -1Q4), (5.5) 
and Eq. (5.3) can be rewritten as 

cos/jQ 4 = aiQi + CX2Q2 + f3 sin fj,Q 3 , (5.6) 
or, renaming (3 as 73 and j3ai by i'ji, we get 

COS//73Q4 = «7iQi + «72<?2 + smfiQ 3 ; = ■ (5.7) 

This equation of course has the required Lorentz covariance if (3.11) transforms as a Lorentz 
vector. Note that without the second step in Eq. (5.2) we would have obtained a different 
equation that does not appear to have the required invariance. 

Let us rewrite Eq. (5.3) as 

e~ iH 1, = [e- if, »Q 3 - iaiQi]^. (5.8) 

On the lattice it becomes 

1>(£, to, t + 1) = e"^ [Fl m ${l + 1, to, t) + F^i - 1, to, t)} + 

+ 73^P^ [tfrnW + 1, m + 1, t) - F} tm W - 1, TO + 1, t)} + (5.9) 
+ l^-^f 1 [ -rfrnW + 1, m ~ 1, t) + Fl m lP(£ 1, TO 1, t)] . 



where the coefficients can be read off from Eqs (4.4)-(4.7): 
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Ft ™ — 



(£ + m+ l)(£-m+ 1) 
(2£ + l)(2£ + 3) 

'(l + m + l)(l + m + 2) 
(2£ + l)(2£ + 3) 

?-m + l)(£-m + 2) 
(2£ + l)(2£ + 3) 



P 2 



P 4 



T7l6 



(£ + m)(£-m) 
(2£- 1)(2£+ 1) ' 

'(l-m)(l-m- 1) 
(2£- 1)(2£+ 1) 

? + m)(£ + m - 1) 
{21- 1)(2£+ 1) 



(5.10) 



By construction, our equation (5.8) produces Eigenvalues with absolute values equal 
to one, and hence we immediately see that the norm of ip is preserved: 



uv(f+i)ii = \\m\\- (5.H) 

This equation can be used as an evolution equation for a quantum state ip. Needless to say, 
our Hamiltonian H has as many positive Eigenvalues as negative ones, and so an eventual 
physically viable theory should have the negative states filled up like a Dirac sea, or else it 
will be impossible to perform ther mo dynamical calculations. As stated before, this would 
go beyond the scope of the present paper. 



6. DISCUSSION 

g 

The introduction of lattice versions of gravity has been proposed before . This par- 
ticuar lattice may also have been considered before. Here we emphasize that the lattice 
structure has been derived from the mathematical expressions for the relation between 
momenta and the Hamiltonian, rather than postulated. Only with the lattice spacing 
precisely defined by the Planck length do the field equations take the form of difference 
equations involving only nearest neighbors on the lattice. 

It should be kept in mind that the lattice structure is more complicated than the 
familiar rectangular or triangular lattices. For a given value of £, the angular momentum 
is limited by the inequality 



\m\ <£. (6.1) 
The distance to the origin, r, is quantized as 
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Concerning the quantum of angular momentum, the reader might have wondered why 
it should be taken as integers (or half-odd-integers) at all. First of all, one could consider 
"anyon" statistics. Now this appears to be difficult in our approach. If energy-momentum 
space is indeed represented by an S2 sphere we will be forced to restrict ourselves to the 
50(3) quantum numbers £ and m, but it is conceivable that more complex theories can 
be constructed; we leave this to other investigators. 

Secondly however, one might argue that we also have to deal with the non-trivial 
holonomies generated by the cusps around the particles. It may seem that angular mo- 
mentum ought to be quantized in units tv/(tv — H) rather than 1. We deliberately postponed 
this complication. The reason is that this topological aspect of space emerges if one particle 
circumnavigates another particle, in other words, it becomes relevant only if we consider 
two or more particles. In the present work only single particle states were considered. 
The center of our coordinate frame is not a particle. It cannot be taken as an infinitely 
heavy particle because such a thing would distort the surrounding space too much. As 
yet we treat the origin as a more abstract reference point in our coordinates. This situa- 
tion is clearly not completely satisfactory. In order to construct a more complete theory 
one will have to addres the many particle situation. ^ One then can restrict oneself to 
describing distances between particles only, while avoiding abstract reference points alto- 
gether. This however will be a major exercise, complicated by the fact that in the polygon 
representation (which we insist on using since it involves proper Cauchy surfaces) shows 
extra contributions to the Hamiltonian from the vertex points between the polygons. We 
suspect that such a treatment will also further clarify the role of the factors cos 11 in our 
expressions (see Eq. 6.1), which presently seem to be rather ugly. 

The equations on our lattice are local in the sense that only nearest neighbors interact, 
and the Cauchy surfaces are well-defined. In contrast, Poincare transformations are in 
general non-local. Related to this is the fact that we have some freedom in choosing the 
details of our lattice. Instead of $2 x S\ we could have picked S3. What happens then is 
explained in the Appendix. The equations also involve nearest neighbors, except that in 
the time direction next-to-nearest neighbors then also appear, and this makes this lattice 
less suitable. The arguments of Sect. 5, in particular Eq. (5.11), do not work in that case. 

The Poincare group in our model is different from the usual one. We still have the 

' An approach towards handling this problem using conformal gauges for space and time is 

9 10 

advocated by Menotti, Seminara and Welling . 
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complete Lorentz group, but the translations are limited to the procedures of adding 
angular momenta using the familiar Clebsch-Gordan coefficients. 

We seem to have hit upon a new terrain of lattice theories where much remains to be 
explored. 
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NOTE ADDED 

Just before this manuscript was mailed a paper by Waelbroeck and Zapata appeared 
in the electronic archives 11 in which also a space-time lattice is discussed in connection 
with 2+1 dimensional gravity. His lattice however is the lattice formed by the imaginary 
parts of the link variables Li as discussed in the Introduction. It is quite distinct from 
what we propose in this paper. 

APPENDIX: THE LATTICE GENERATED BY S 3 . 

The algebra of the spherical harmonics on £3 is generated by first considering rotation 
operators L a and M a , a = 1, 2, 3, in a four-dimensional space with coordinates (qi, ... , q 4 ): 

L a = ~ie ahc q b d c ; M a = i(q 4 d a - q a d 4 ) . (Al) 
They obey the commutation rules 



[L a ,L b ] = ie abc L c ; (A.2a) 
[L a ,M b ] =ie abc M c] (A 2b) 

[M a ,M b ] =ie abc L c . (A2c) 



Construct the self-dual and the anti-self-dual combinations, 



L L a =\(L a + M a ) ; L? = \(L a - M a ) . (A3) 
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We have 



[L£,L fe L ] = ie ahc L L c ; [L? , L? ] = ie a6c Lf ; [Lj,L«] = 0. (A4) 
Since L • M = we find 

L l2 = L r2 =i(L 2 + M 2 ) = 

The simplest non-trivial function is fn(q) = q^, which has 

dlU = ; (qMU = f, 

so that we read off from Eq. (A. 5) that it has t=\. Its four components are represented 
by rriL = ±|, rrtR = These functions / M are then to be considered on the sphere 

\q\ = 1. By multiplying these together we can construct functions with all other (integral 
or half integral) £ values, whereas on the other hand it is easy to convince oneself that 
every set of numbers (£, rriL, tor) with |m^| < £ and \mn\ < £ corresponds to one and only 
one spherical harmonic function (up to a phase factor). These numbers £, tul and tur 
therefore form a complete set of coordinates. 

We have a freedom in our choice for orienting the coordinates in the g M space. 
How one does this turns out to make little difference; in all cases the £, ttil, ttir do not 
directly correspond to all space-time coordinates. A convenient choice is: 

Qi = qi , Q2 = qi , Q3 = 94 , Qi = qs ■ (A.7) 

Then the variable conjugated to the Hamiltonian (3.10) is 

t = -M 3 = m R -m L . (A8) 
and from Eqs. (3.2) and (A.l) we have 

x cos fx = L 2 ; y cos \i = —Li, (^.-9) 

Consider an arbitrary function of space-time, i.e., one that does not obey any equation 
of motion such as a mass shell condition. For such a function the three coordinates x, y, 
and t can all three be regarded as operators. They are quantized according to 



£{£+!) 



(A.5) 



(A.6) 
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%i = n>i I cos fj, ; t = 713/ cos \x , 



(A10) 



however, these three numbers n\ are non-commuting. We have 

[x,y] = i(cos 2 /i)L 3 , 
[y,t] = -i(cos /z)M 2 , (All) 
[x,t] = —i(cos /j,)Mi . 

The operators £1,2,3 are the usual generators of Euclidean rotations among the Mj. Thus 
only the operators t, L 3 , and £ are diagonalized simultaneously. I roughly corresponds to 
'total Euclidean distance from the origin'. If we would limit ourselves to integral I we find 
that at even time — mi is even, therefore + is even. This would imply that 
at even time parity of all functions f(x, y) is even: f(x, y) = f(—x, —y), and at odd times 
parity is odd. If one wants to avoid that, half-odd-integer values of I should be admitted 
together with the integral ones. 

The Li represent Euclidean rotations in space-time. Lorentz invariance is not a prop- 
erty of the coordinates but of the mass shell condition (3.3). Lorentz transformations are 
defined in terms of Eq. (3.11) in energy-momentum space. They map the sphere onto itself, 
after which decomposition in spherical harmonics will give the new coordinates. Since the 
L's and M's of Eqs (A. 2) are the generators of Euclidean rotations they do not transform 
into combinations of each other under Lorentz transformations. Therefore the Lorentz 
group will act in a fairly complicated way on our coordinates; this statement holds both 
for the £3 case and for the S2 x Si case treated in the text. 




Fig. 3. The mass shell: a) in ordinary 2+1 dimensional particle theory, b) in 2+1 di- 
mensional gravity for the £3 case. Energy and momentum live on an S3 sphere. In 
both cases the two space axes are drawn as one. 



19 



Now consider the mass shell condition, Eq. (3.3), which in the S 2 x Si case had to be 
rewritten as Eq. (4.2). In order to be able to process it in £3 we have to rewrite it as 



Ql cos 2 /! = Q\ + Q\ + Ql sm 2 /j (A12) 

(compare Eq. (5.2)). In ordinary flat space it is replaced by Eq. (3.3a), which corresponds to 
the Klein-Gordon equation in coordinate space. In Fig. 3 we show how the mass shell looks 
on the S3 sphere. What is the equation generated by (A. 12) in our space-time lattice? The 
functions Q^Q V form a representation of SO {A) obtained from the symmetrized product 
of two £ = j representations Q M : 



4x5/2 = 10 = 9 + 1, (A13) 

where the 1 is of course the scalar Q 2 , which is fixed to be 1 on the sphere. The 9 is formed 
from the states £ = 1, |mz,| < 1, \mp\ < 1. Multiplying a function at the point £,m,L, tur 
with this representation leads one to a neighboring point I + A£, mi + Am^, rriR + Am^, 
with \ A£\ < 1, Am_L| < 1 and |Am^| < 1. We see that, as in Sect. 4, the Klein-Gordon 
equation is replaced by a difference equation on the lattice. We also see that only points 
differing by integral amounts of £, ttil and rtiR are connected, so that the space-time lattice 
breaks up into two mutually non-interacting systems. This corresponds to the conservation 
of parity. 

Let us consider the action of the operators (or equivalently the q^) in somewhat 
more detail. From Eqs. (A.l) and (A. 3) we have 



L 3 h = ^i{q 2 di -qid 2 ± <?4<9 3 =F <?3<9 4 ) . O 4 - 14 ) 

In order to expand a "state" ^ a^q^ in terms of eigenstates ||, m^, to#), we solve the 
equations 



1 

L 3 R (a ■ q) = ~i(q 2 ai - <?i«2 ± ?4«3 =F <?3«4) 

(A15) 

= m l (otiqi + a 2 q 2 + a 3 q 3 + 014(74) , 

Ft. 

Indicating the states 1 1, ra^ iur) (where of course all m values are ±|) by the short-hand 
notation an), where the cr's denote the signs of the m's, we find: 
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v^| + +> 
>/2| + -> 
>/2|-+) 
v^|--> 





«2 


«3 


CK4 






1 


% 








\/2gi = 


I + - I"") 








-1 


—i 


V2q 2 = 


-i| + +) - i\ - - 








-1 


% 


v^ttj = 




-1 


z 








V2q* = 


i\ + -> - i\ - +> 



(A16) 



The relative signs here have been fixed by considering the (left and right) angular momen- 
tum raising and lowering operators. 

Now that we know their left and right quantum numbers we are in a position to 
calculate the action of the operators on the states \£,m,L, thr). The m quantum numbers 
just add up, whereas the total angular momentum numbers £ change by one half unit. Thus 
we can state that 



\<r L ,<TR)\e,m L ,m R ) =cC t ''mjH 5, m L + \<r L , m R + \(Jr) + 

where a and (3 are coefficients that can now be calculated. The technique to be used is a 
standard manipulation with Clebsch-Gordan coefficients, similar to what we did in Sect. 
4, but we will not repeat it here. 

We see confirmed that the equation (A. 12), in which each term corresponds to acting 
twice with operations as described in Eq. (A. 17), amounts to a difference equation on the 
(£, rriL, trr) lattice. In general this equation takes the form: 



^2 F(£,m L ,m R ,a e ,a L ,aR)ip(£ + ae,m L + a L ,mR + aR) = 0, (A. 18) 

where all variables (Je,L,R take the values ±1 only. Now, in view of Eq. (A. 8), we see that 
to describe the complete space of solutions to this equation one needs the Cauchy data on 
several consecutive "Cauchy planes": t = t\ + 1, ti, . . ., t\ — 2, after which the values at 
t = tt+2 follow. 

The Dirac equation again takes the form 



hiQi + 73<24Cos / u + (23sin/i = 0, (A. 19) 

where 71,2,3 can be taken to be the three Pauli matrices. Multiplying the equation with 
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HiQi + H2Q2 + 73Q4 cos \i - Q 3 sin \x , 



(A.20) 



reproduces (A. 12). Note that Lorentz covariance follows from the usual arguments, 
using the transformation law (3.11). 

Consider a solution of the difference equation (A. 19). In the notation of Eq. (A. 16) 
our equation (A. 19) takes the form 



(72 + *7i)l + +) + l72-*7i)|- -)+^73e l73M | + -)-i73e- l73At |-+) = 0. (A21) 

Since t = uir — rriL (Eq. A. 8), the last term advances time by one unit and the third term 
retards by one unit. The first two terms stay on the same time sheet. One can never have 
a positive and conserved probability density for such a wave equation. 

In conclusion, besides giving mathematically more cumbersome difference equations, 
the S3 lattice is also physically less appealing than S2 x Si. Ultimately, however, both phys- 
ical systems are mathematically equivalent, since in momentum space the field equations 
are identical. 
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